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Above 54.4 eV, two-photon double ionization of helium is dominated by a sequential absorption 
process, producing characteristic behavior in the single and triple differential cross sections. We 
show that the signature of this process is visible in the nuclear recoil cross section, integrated over 
all energy sharings of the ejected electrons, even below the threshold for the sequential process. 
Since nuclear recoil momentum imaging does not require coincident photoelectron measurement, 
the predicted images present a viable target for future experiments with new short-pulse VUV and 
soft X-ray sources. 
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Kinematically complete experiments on one-photon 
double ionization of the simplest atomic jl] and molec- 
ular [5] targets, coupled with state-of-the-art non- 
perturbative theoretical studies, have yielded fundamen- 
tal information and insight into the nature of electron- 
electron correlation. With the advent of new classes of 
short-pulse, high intensity free-electron laser [3l and high- 
harmonic generation ^ light sources that operate in the 
VUV and soft X-ray regimes, one looks forward to sim- 
ilar experiments involving the nonlinear process of few- 
photon multiple ionization [51 With plans for such 
experiments currently underway, it is important to un- 
derstand what kinds of new phenomena can be studied 
and the measurements that are most likely to be success- 
ful in revealing new effects "1_ ,. 

Helium offers an interesting case in point. The energy 
required to doubly ionize helium is 79.0 eV, the sum of 
the first (24.6 eV) and second (54.4 eV) ionization ener- 
gies. Double ionization of helium by two-photon absorp- 
tion therefore requires a minimum photon energy of 39.5 
eV. For photon energies between 39.5 and 54.4 eV the 
process is necessarily direct or nonsequential (NSI), ie. it 
requires simultaneous absorption of two photons to strip 
two electrons from the atom and the process is expected 
to be sensitive to electron-electron correlation. For pho- 
ton energies above 54.4 eV, sequential ionization (NSI) 
is possible: He"*" is produced by absorption of one pho- 
ton, followed by absorption of a second photon to produce 
He"'""'" . Although NSI can compete with SI at these higher 
energies, SI, which is an essentially uncorrelated process, 
will dominate. In NSI, the excess energy {2hhj — 79 eV) 
can be shared continuously between the photoejected 
electrons, whereas with SI, we expect to see photoelec- 
tron energies sharply peaked around — 24.6 eV and 
hw — 54.4 eV. Experimental studies E], and indeed 
most theoretical treatments, of this process have to date 
focused on energy and angle-integrated quantities. In 
contrast to the one-photon, double ionization case, kine- 



matically complete experiments that could distinguish 
between NSI and SI have yet to appear. Our purpose 
here is to demonstrate, using the results of precise quan- 
tum mechanical calculations, that since the angular dis- 
tributions produced in two-photon NSI and SI are very 
different, they leave a clear signature in the differential 
nuclear recoil distributions, even when integrated over all 
possible energy sharings of the ejected electrons. We will 
show that the signature of SI can be seen in the pho- 
toelectron energy-integrated nuclear recoil cross sections 
at energies several eV below the 54.4 eV threshold where 
sequential ionization is still a virtual process. Since these 
distributions do not require direct detection of the photo- 
emitted electrons, they present a attractive target for an 
experiment that could demonstrate a clear signature of 
sequential versus nonsequential double ionization. 

The triple differential cross section (TDCS) for two- 
photon double ionization is defined as 

da 27r(27ra)2 2 /-.n 
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where ki and k2 are the momenta of the photoelectrons, 
Lj is the photon frequency, m is the electron mass, a is the 
fine-structure constant, and Ex = k'l/2 is the energy of 
one of the electrons and thus defines the energy sharing. 
The ionization amplitude /(ki,k2,a;) is in turn given by 

/(ki, k2, Lo) = (vl/^^ \^,[Eo + hLO-H + ze]- Vl'&o) . (2) 

H is the atomic Hamiltonian, $0 is the initial state of 
the atom with corresponding energy Eq, and "i!^^ is 
the full momentum-normalized scattering wave function, 
with incoming boundary conditions corresponding to two 
free electrons,. For polarization e, the dipole operator in 
the velocity form, /i, is defined in terms of the momentum 
operators, Pi, for the two electrons by /i = e • pi -I- e • p2. 

The accompanying momentum recoil, Q, of the nucleus 
due to the ejection of two electrons of momenta ki and 
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k2 is 



Q = -(ki+k2) 



(3) 



At a given photon energy, we can therefore define a nu- 
clear recoil cross section, differential in energy sharing 
and the angular dependence of Q, by 



da 



d^QdEi 



^ dfli da 



da 



daidQ.2dE] 



■(53(Q + ki+k2) 

The integral of this quantity over energy sharing, 
da/d^Cl — J dEi da /d'^QdEi, carries a signature of the 
sequential ionization process even below its threshold, as 
we will see shortly. 

We have recently shown lOJ that task of calculating 
the two-photon, double ionization amplitude, which re- 
quires in principle the exact wave function for three-body 
Coulomb breakup, can be simplified by using the method 
of exterior complex scaling (ECS) which avoids the ex- 
plicit imposition of asymptotic three-body boundary con- 
ditions [11] . We begin with the coupled driven equations 
in the Dalgarno-Lewis form of second order perturbation 
theory [12 that describe the absorption of two photons 
by a system initially in state <i>Oi 



(So + ^-i^)*r(ri,r2) = M$o 



(^0 + 2hiu - i/)^'r(ri,r2) = fi^ 
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which must be solved with pure outgoing wave bound- 
ary conditions for the wave functions "iff^ and ^'|'^ . With 
ECS, Eqs. ^ and ^ are solved numerically on a dis- 
cretized grid in a large but finite region of coordinate 
space extending to some Rq, where the outer boundary 
conditions are obviated by rotating the radial coordi- 
nates beyond that point into the complex plane. Having 
solved the driven equations, the ionization amplitude can 
then be extracted using a surface integral that involves 
a pair of testing functions (r) which are momentum- 
normalized, one-electron Coulomb functions with nuclear 
charge Z^2, in the case of helium [T^ : 

/(ki,k2,w) = i I (^k;(ri)Vk;(r2)VvI/-(ri,r2) 

-vI/-(ri,r2)V^k:(ri)V'k;('^2)) -dS 

(7) 

The integral is evaluated over a finite hypersphere whose 
radius is less than Rq. No approximation concerning the 
final state has been made in this formalism because the 
testing functions, ip^iv) , merely extract the double ion- 
ization amplitude from the final outgoing wavefunction 
^'1'^, and electron correlation is treated completely in ^'l'^ 
as well as in the initial state $0 in this approach. 

But there is a technical difficulty must be addressed. 
For photon energies greater than the first ionization po- 
tential of the atom, ^f^'^, the solution of Eq. (|5|, will have 



undamped outgoing wave behavior on the real portion of 
the grid. Since ^'f'^ serves as the source term in Eq. (|6|, 
^'l'^ will not converge with increasing Rq. As explained 
if ref. |10j , we can solve this problem by adding a small 
positive imaginary part to to in Eq. ^ only, rendering 
"iff^ square-integrable, and then numerically extrapolat- 
ing the results to real w. 

'i'f^ and ^'|'^ were expanded in a product basis of spher- 
ical harmonics, giving a set of coupled two-dimensional 
radial equations that were discretized using a finite el- 
ement, discrete variable representation (DVR) [Hj and 
solved on parallel computers using sparse matrix meth- 
ods. For the results reported here, we used partial waves 
up to Z = 9 and radial grids with real parts extending to 
160 bohr on a side, with finite element boundaries start- 
ing at 5 bohr and then spaced 10 bohr apart, discretized 
with 18th order DVR. Calculations were performed for a 
range of complex lo values between Im(a;) — 0.500 and 
Im(a;) = 0.05 hartrees and the individual partial wave 
amplitudes were then extrapolated to real uj. 

Figure [T] shows the single differential cross section 
(SDCS), which is obtained by integrating the TDCS de- 
fined in Eq. ([ij over the angles fii and 0,2, at three differ- 
ent photon energies. The SDCS is seen to be a relatively 
flat function of energy sharing below 50 eV. The signa- 
ture of SI becomes evident above 50 eV where the SDCS 
begins to rise at the extremes of energy sharing until the 
sequential limit is reached at 54.4 eV. Above that en- 
ergy, it has well defined peaks at Ei — huj—bAA eV and 
El = eV. Our calculated SDCS at 58 eV has 

finite peak heights because it was extrapolated from cal- 
culations carried out on a finite grid. These peaks would 
become singularities in the limit of an infinite grid, and 
are a fundamental feature of the cross sections derived in 
lowest order perturbation theory. The apparent widths 
of the peaks, however, are not a consequence of the finite 
grid, as explained below. 

One also sees a striking change in the TDCS as we 
move from regions dominated by NSI processes to those 
where SI dominates. Figure [2] shows contour plots of the 
TDCS, as a function of energy sharing and scattering an- 
gle for one of the ejected electrons, at photon enertgies 
of 44, 52 and 58 eV. The scattering angle of the sec- 
ond electron is fixed at 30°. At 44 eV, the TDCS has 
its maximum value when the two electron are ejected 
back-to-back, independent of the energy sharing, which 
is expected for a correlated, NSI process. At 52 eV, the 
angular dependence is more complicated, but one still 
sees a propensity for back-to-back ejection, as well as a 
peak near 180° for extreme unequal energy sharing. At 
58 eV, the TDCS is strongly peaked near and 180° at 
the energy sharings corresponding to SI and is uniformly 
small elsewhere. 

As we showed earlier |10j, we can get a qualitatively 
correct description of sequential ionization with a simple 
analytic model that ignores correlation in the initial, in- 
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FIG. 1: Single difTerential cross sections at 44, 52 and 58 eV. 




FIG. 2: Triple differential cross sections, as a function of en- 
ergy sharing and angle, at (top to bottom) 44, 52 and 58 eV, 
showing the signatures of nonsequential and sequential ion- 
ization (see text). The direction of one electron is fixed at 
30° to the polarization direction while the other is varied. 



termediate and final states. That model, whose details 
are derived elsewhere |15j . gives the following expression 
for the TDCS, 



h 
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dEidnidQ.2 
Eq + huj — eis — El 



cos2(6ii)cos2(6'2) 

^a^'^+{Ei)a^-iE,) 
Eo + hiu ~ eis - E2 



(8) 



where a^'^ and 0-^°^ are the photoionization cross sec- 
tions of helium and He+, respectively, and eis is the en- 
ergy of He"*". When integrated over fii and ^2, Eq.Q 




FIG. 3: Nuclear recoil cross sections, integrated over all pho- 
toelectron energy sharings, at 44, 52 and 58 eV. The polar- 
ization vector is chosen to lie along the z-axis. 



gives an SDCS proportional to the square of the term 
in square brackets. Thus the model shows that when the 
photon energy exceeds 54.4 eV (eis) the TDCS and SDCS 
diverge at the ejected electron energies corresponding to 
SI, the singularities separated by the difference between 
the ionization energies of He and Hc+. The model also 
explains the rapid rise of the SDCS just below the SI 
threshold at extreme unequal sharing, since an appar- 
ent width (determined by the size of the numerators in 
Eq.(|8])) is associated with the singularities above the SI 
threshold which causes the SDCS to rise above the NSI 
background . As mentioned above, these singularities 
are consequences of treating sequential ionization in low- 
est order perturbation theory, and not of the simplifying 
assumptions used in deriving the model. 

The simple model also shows that the angular depen- 
dence of the TDCS is the product of two uncorrelated 
dipole distributions for each ejected electron. This be- 
havior is clearly seen in the calculated TDCS at 58 eV 
which shows peaks at O2 — 0° and 180°. The calculated 
TDCS can be used to compute the momentum recoil im- 
parted to the nucleus. 

Figure |3] shows the calculated nuclear recoil cross sec- 
tions, integrated over all energy sharings, at the three 
photon energies previously considered. The integration 
m Eq.(|4| requires some effort because, for example, there 
are ranges directions of ki, for which no k2 exists that 
satisfies Eq.([3|. Details about the numerical evaluation 
of da / d^CldEi will be given elsewhere |15j. The gen- 
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eral shapes of these cross sections can be understood on 
the basis of the previously described TDCS. For NSI, 
the two electrons are preferentially ejected back-to-back, 
which means comparatively little momentum is imparted 
to the nucleus. But the imparted momentum for back-to- 
back ejection is exactly zero only for equal energy shar- 
ing and we are integrating over all energy sharings, so 
we would expect a roughly isotropic distribution peaked 
about Q = for a pure NSI process. For SI, however, the 
electrons are uncorrelated and, because of their cos^{9) 
angular dependence, are preferentially ejected along the 
polarization axis, i.e. at 0° or 180°. We therefore expect, 
for a pure SI process, that the distribution of nuclear re- 
coil momenta will show four peaks along the polarization 
axis at IQI = ±(|ki| + jkal) and |Q| = ±(|ki| - jkal). 
These general features are indeed seen in the calculated 
cross sections. The nuclear recoil cross sections at 58 eV 
show four clearly defined rings (whose peak heights we 
must emphasize would diverge if we could use an infi- 
nite grid). Even at 52 eV, we can see a clear signature 
of SI in the nuclear recoil cross section where there are 
two distinct rings at the extremes of momentum trans- 
fer along with two secondary wings developing inside the 
prominent wings - the signature of "virtual SI" . 

In summary, we have shown that sequential ionization 
leaves a clear signature in the nuclear recoil cross sec- 
tion, even when it is integrated over all energy sharings 
of the ejected electrons. This signature is clearly visi- 
ble at 52 eV, which is 2.4 eV below the SI threshold. 
The nuclear recoil cross section is found to be a very 
sensitive indicator of SI. We see "virtual SI" structure 
in the nuclear recoil cross section even at 44 eV, while 
the SDCS at the same photon energy is flat. The energy- 
integrated nuclear recoil can be measured in a double ion- 
ization experiment using, for example, the COLTRIMS 
technique [16]. Such observations do not require detec- 
tion of the electrons at all and thus avoid completely the 
need for coincidence measurements. 
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